


 % Quit rate
 %%%%%%%%%%%%
 % Natural wastage
quit_nw = s*lambda*ones(obs_nw,1);
 %
 % Full replacement
quit_fr  = s*lambda ./ (1 + aux*s*lambda*log(grid_fr/mM));
 %
 % Net Expansion
delta_mR = s*lambda  / (1 + aux*s*lambda*log(mR/mM));
auxconst = ( ((1-omega1)/(alpha*c))*mR - (r+delta_mR + ((omega0+r*C)/c)) )*(mR^(-1/(1-alpha)));
quit_ne  = ((1-omega1)/(alpha*c))*grid_ne - (r+((omega0+r*C)/c)) ...
           - auxconst*(grid_ne.^(1/(1-alpha)));


 % Hiring rate in Net Expansion region
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
qq_ne         = (ushr/(1-ushr)) + G_ne;
deltaprime_ne = ((1-omega1)/(alpha*c)) - auxconst*(1/(1-alpha))*(grid_ne.^(alpha/(1-alpha)));
eta_ne_x_g_ne = -deltaprime_ne.*qq_ne;
eta_ne   = eta_ne_x_g_ne./g_ne;
